Dynamical Phases of Driven Vortices Interacting with Periodic Pinning 
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The finite temperature dynamical phases of vortices in films driven by a uniform force and 
interacting with the periodic pinning potential of a square lattice of columnar defects are investigated 
by Langevin dynamics simulations of a London model. Vortices driven along the [0,1] direction and 
at densities for which there are more vortices than columnar defects (B > B^) are considered. At 
low temperatures, two new dynamical phases, elastic flow and plastic flow, and a sharp transition 
between them are identified and characterized according to the behavior of the vortex spatial order, 
velocity distribution and frequency-dependent velocity correlations. 
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There is a great deal of current interest in the study 
of dynamical phases of vortices driven by an external 
force and interacting with various arrays of pinning cen- 
ters. One possibility is periodic pinning resulting from 
a lattice of artificial defects. Understanding the dynam- 
ical phases in this case is of interest because it may be 
possible to observe them in superconducting films with 
periodic arrays of holes, magnetic dots and columnar de- 
fects (CD). Several techniques have been developed to 
fabricated these films and studies of vortex dynamics in 
them have been reported Jl],^[ . Theoretical workers have 
carried out investigations, mostly numerical, of driven 
vortices under periodic pinning However, as dis- 

cussed here, several questions remain open. 

To be specific, this paper considers two-dimensional 
vortices interacting with a columnar defect lattice (CDL). 
The motion of N v such vortices is assumed to be governed 
by Langevin equations for massless particles, which for 
the l-th vortex reads [Ml , 



dri 
77— = F 
' dt 



-cd 



(1) 



where r\ is the friction coefficient, 



— ^2j^ l=1 ViU vv (yi — Tj) is the force of interaction with 

other vortices, F^~ cd = -^Zn^i uv ~ cd ( T i ~ R ) is the 
force of interaction with the CDL, R denotes the CDL 
positions, is the driving force, and Ti is the random 
force appropriate for temperature T. 

In the absence of driving (F^ = 0), the equilibrium 
phases of this system ( hereafter called zero-drive state) 
are studied (at low T) in Refs. §-0. 

For F^ 7^ the simplest situation occurs at very 
high drives (jnj. In this case the vortices center of 
mass (CM) move with velocity = F^/rj and in the 
CM frame of reference, defined by rj = r$ — V^t, the 
equations of motion are as in Eqs. (Q), but without 
the driving term F^ and with the vortex-CDL force re- 



E Q (-^Q)^" C(i (Q)e iQ - r; e iQ - v ^, where Q denotes the 
CDL reciprocal lattice vectors and U v ~ cd (Q) is the 
Fourier transform of U v ~ cd (r). For high enough drives 
the Fourier components in F^ v ~ cd \t) for which Q • ^ 
oscillate fast, thus having a negligible effect on the 
vortex trajectory [jn]. The vortex-CDL force in the 
CM frame reduces then to the static one obtained by 



summing the Fourier components in F 



t{v-cd) 



(t) with 



placed by the time-dependent interaction F 



f(v—cd) 



(*) 



Q • Vd = (or Q _L V^). This force is just that result- 
ing from the potential obtained by averaging JJ v ~ cd (Y) 
in the direction of drive. Thermal fluctuations lead, af- 
ter sufficient time, to relaxation into the thermodynamic 
equilibrium state of the N v vortices interacting between 
themselves and with this average potential. This state is 
referred to here as the infinite- drive state. 

Since at low T the infinite-drive and the zero-drive 
states are generally different, reordering of the vortices 
must take place as the driving force is varied. This pre- 
dicts the existence of at least two dynamical phases and, 
possibly, of one dynamical phase transition. One novel 
result obtained in this paper is to show that this pre- 
diction is correct. Studies of vortices interacting with a 
CDL carried out at T = find a rich variety of phases, 
but no reordering to the infinite drive-state ^,||. 

For driven vortices interacting with a random pinning 
potential, the existence of reordering and of a dynamical 
transition as the driving force is varied are well estab- 
lished ]l2|. In this case the pinning potential averaged 
in the direction of drive remains random in the direction 
perpendicular to it, and has non-trivial effects on the dy- 
namical phase diagram, as extensively discussed in Refs. 
|Jl3f . For motion on a periodic pinning potential, numer- 
ical simulations reported in Refs. Jl4|] show that vortex 
reordering to the infinite- drive state only occurs if suf- 
ficient thermal fluctuations are present. Otherwise the 
vortices are trapped in a metastable state. 
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The numerical simulations reported here investigate 
the properties of the moving vortices steady-states as a 
function of the driving force, for a range of T beginning 
at a value well below the infinite-drive state melting tem- 
perature, T m , and extending up to T m . The simulations 
are initialized with the vortices in the infinite-drive state 
and large enough F d . In subsequent runs, at the same 
T, F d is progressively decreased. The main conclusions 
reached by this approach are that for T < T m a dynami- 
cal phase with the spatial symmetry of the of the infinite- 
drive state, elastic flow and long-range time order exists 
for Fd > F C (T). At Fd = F C (T) a transition separates it 
from another dynamical phase with distinct spatial order 
and plastic flow. It is found that the transition is sharp 
at low T and that F C (T) increases as T — > T m . 

The details of the model are as follows |H| ]. Vor- 
tices and CD are placed on a square lattice subjected 
to periodic boundary conditions (hereafter called the 
space lattice), with N square primitive cell of dimen- 
sions d x d, oriented with the sides in the x and y- 
directions. The vortex has a core of linear dimension 
d v > d (d v ~ 2£(0)). The interactions between vortices 
and between vortices and CD are chosen to model films in 
the London limit. The vortex- vortex interaction poten- 
tial is a screened Coulomb one appropriate for the 
space lattice, with a short distance cutoff to account for 
the vortex core [Q. That is, U vv (r) is the lattice Fourier 
transform of U vv (k) = 4tt 2 Jexp (- n 2 / n 2 ) / (n 2 + A~ 2 ), 
where J = (cj^dy / ?>2tt 2 \ 2 ) is the energy scale for vortex- 
vortex interactions |l5| , A is the penetration depth, 
k 2 = 4 sin 2 (k x d/2) +4 sin 2 (k y d/2), k c = 2 sin (ird/ 2 d v ) 
is the cutoff in k-space, and A is the screening length 
(A > A). The CD-lattice has N c d sites arranged on 
a square lattice, with lattice constant a cd , and is com- 
mensurate with the space lattice. The vortex-CD po- 
tential is chosen with depth U v ~ cd (0), range R c d > d v 
and a spatial dependence that gives square equipoten- 
tials and an attractive pinning force of constant mod- 
ulus F p =| U v ~ cd (0) | /R cd . That is, U v - cd (r) = 
U v - cd (0)(l- | x | /R cd ) for | x \< R cd: - | x \< y <| x |, 
and U v - Cd (r) = U v - cd (0)(l- | y | / R cd ) for | y \< R cd , 
- | y |< x <| y |. 

The algorithm for the simulation of this model follows 
the usual procedure |T(j|| . In the results reported next 
space lattices with N = 256 x 256 and 512 x 512 sites 
are used to accommodate a CDL with N cd = 64 and 
256, respectively, both with a c d = 32d. Other parame- 
ters are fixed at d v = 4d, A = 160d, U v ~ cd (0) = -J, 
Red — 3d v = 12d, from which it follows that F p = J/12d. 
The vortex systems studied have N v = 1.2hN cd and 2N c d 
or, in terms of the magnetic induction, B = 1.25P^ and 
B = 2P^, where B^ = </>o/a 2 d is the matching field. Typ- 
ical run times are 2 — 5 x 10 6 steps of 10 _2 r (r = r]d 2 / J is 
the unit of time), after the steady state is reached. The 
reported results represent the average over several differ- 
ent realizations of the random force(~ 10). The driving 
force is assumed to be along y (CDL (0, Indirection). 
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FIG. 1. a) Structure function for B — 2B$ and 
T/T m = 0.1. Squares (circles): k= typical infinite (zero)- 
drive state reciprocal-lattice vector of smallest modulus. In- 
set F c vs. T. b) Corresponding P(r) for: 1) infinite- drive 
state, 2) elastic flow phase at Fd = 0.7 F p , 3) plastic flow 
phase at Fd = 0.6P P , 4) zero-drive state for Fd < 0.45P P . 
The origin r = is at each picture center. Gray scale in each 
picture is proportional to P(r), but scales in different pictures 
are unrelated. 

To characterize the dynamical phases the following 
quantities are calculated in the steady-state regime: i) 
Time average of the individual vortex velocities, Vj (j = 
1,...,N V ), their mean or the center of mass (CM) ve- 
locity, V cm = "ZZjVj/Ny, and root-mean square devia- 
tion AV^ m = yjEj(vja-Vc™) 2 /Ny, where a = x,y. 

ii) Time-averaged density-density correlation function, 
P(r), and its Fourier transform, the structure function, 
S(k). hi) Time-dependent correlation functions for the 
CM velocity, C^ m (t) = (V^ m (t + s)V^ cm (s)) s , and for the 
vortex velocity autocorrelation function averaged over all 
vortices, C S J (t) = ^- ^Z i (v i a(t^s)v ia (s)) s , where () s de- 
notes average with respect to time s, and the respective 
Fourier transforms C^ m (a;) and C^/ (u). Physically, P(r) 
is proportional to the probability of finding a pair of vor- 
tices separated by r, whereas C^ m (uj) is proportional to 
the noise power spectrum. 

The possible infinite-drive states and the correspond- 
ing T m are determined by equilibrium Monte Carlo sim- 
ulations of the model described above, with the vortex- 
CDL interaction potential replaced by its average in the 
^-direction. This is a periodic washboard potential, con- 
sisting of grooves of width 2R c d running along the y- 
direction and separated by a cd along the x-direction. In- 
side the groove the vortex is attracted to its center by 
a force of constant modulus F < F p in the x-direction. 
For B > B§ these states have only two possible spa- 
tial orders: incommensurate (essentially triangular) at B 
sufficiently large, and commensurate with the ^-direction 
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periodicity at lower B. The precise B value at which the 
ground-state changes from one type to the other depends 
on a non-trivial way on the model parameters. It is found 
that, for the above described model parameters, two typ- 
ical such infinite-drive states are the vortex-lattices (VL) 
for B = 2B0, with the -P(r) shown in Fig. 0.b.l, and for 
B = 1.25^0 with the vortex arrangement shown in Fig. 
||.a. The latter can also be seen as consisting of pairs of 
vortex chains located within the grooves and displaced 
relative to one another by half intra chain spacing. It is 
found that the VL for B = 2B§ melts to a vortex liquid 
at T m = 0.9J, and that the VL for B = 1.25P first 
melts to a smectic phase at T/T sm ~ 0.85 J and then 
to a liquid at much higher T. The results reported here 
for B = 1.25^0 are restricted to T/T sm ~ 0.1. Monte 
Carlo simulations are also used to determine the zero- 
drive states for these B values. For B = 2P^ a VL com- 
mensurate with the CDL and having one extra vortex 
at the center of the CDL primitive unit cell is obtained, 
in agreement with previous results PJl0|]. This state's 
P(r) is shown in Fig. [l].b.4. For B — 1.251^ a complex 
commensurate VL with several vortices in the unit cell is 
found. 
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FIG. 2. B — 1.25^0 : a) Full circles: vortex positions 
for infinite- drive state. Full line: VL primitive unit cell. 
Open squares: CDL. b) CM velocity Vy 771 (units d/r) and 

5 Cm = AV y C7n /Vy Cm . 

The results of the dynamical simulations carried out 
here, to be discussed in detail next, identify two dynam- 
ical phases, elastic flow and plastic flow, and a transition 
between them around Fd = F C (T). It is found that for 
B/B^ = 1.25 and T/T m ~ 0.2, F c /F p ~ 1.25 and that 
for B /Bff, = 2 the F C (T) curve is that shown in Fig. [lja. 
In both flow regimes the time-averaged velocity of all vor- 
tices is in the direction of drive. The Vy m (Fd) curves ( 
V-I curve ) at low T are depicted in Figs, g.b, ||.a. The 
detailed properties of these phases are as follows. 

i) Elastic flow: F d > F c , low T, both B. The mov- 
ing vortices remain localized with respect to each other 
around the infinite-drive state relative equilibrium posi- 
tions. This is suggested by P(r) shown in Figs, [jj.b.l, 
|l|.b.2 and [I|a. These consist of isolated spots with the 
symmetry of this state. As Fd — > F c these spots be- 
come larger, indicating increase in relative vortex mo- 
tion in the CM frame. The time-averaged velocity of 



AV^/VJT 1 (a = x,y) is small for F d > F c ( Figs. 
||.b, and[|.a). Similar results are found for AV^ 771 /Vy 771 . 
Vortex motion consists of a translation in the direction 
of drive with Vy 771 and of periodic oscillations with fre- 
quency uo crn = 27rVy m /a c d and higher harmonics around 
the infinite-drive state equilibrium positions. These os- 
cillations are suge; ested by Q m (cj) and C S J (uo) which are 
found to have sharp peaks at uo cm and peaks with smaller 
amplitudes at some harmonics as illustrated in Fig. ||.b. 
This frequency coincides with the oscillation frequency 
of the vortex- CDL interaction in the frame moving with 
the CM. This suggests that the periodic vortex motion is 
the (non-linear) elastic response of the VL to this interac- 
tion. The frequency dependence of C^ m (c^) and C S J (cu) 
suggest long-range time order. 

ii) Transition region: Fd ~ F c , low T, both B. Vor- 
tex spatial order changes in a small Fd interval. This 
is seen both in P(r) and in 5(k), for k equal to high- 
drive ground-state reciprocal lattice vectors, (Fig. [I]). 
The value of F c is estimated as the mid point of this 
interval, as indicated in Fig. |l|. This rapid change in 
the vortex spatial order suggests a discontinuous jump in 
5(k). However, the present data cannot rule out a sharp 
crossover. In the same Fd interval, the Vy 771 x Fd curves 
changes slope, as seen in (Figs. [|b, |[a) and the sharp 
peaks in C^ 771 ^) and C^/ (cj) are found to disappear. 
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5 C m = AV y crn /V y crn . b) Correlation function (u) at F d /F p 
values indicated in the figure. The ordinate axis (arbitrary 
units) in each curve is scaled to show all curves in the same 
graph, uo denotes uj crn at Fd/F p = 1.0. Peaks are centered 
at co cm in all curves. 



iii)P/as^c flow: F d < F c , low T, B = 2P^. Plastic 
flow sets in for Fd just below F c and relative motion in- 
creases as Fd decreases, as evidenced by the sharp change 
in slope and continuous growth in AVy 771 /Vy 771 shown in 
Fig. ||.a. The vortices become pinned in the zero-drive 
state for Fd < 0.45P p . It is found that all vortices become 
pinned simultaneously at Fd = 0.45P p . The vortices are 
disordered for Fd just below F c (Fig. |l|.b.3) and, as Fd 
further decreases, the vortices show increasing order in 
the zero-drive state symmetry, as suggested by the S(k) 
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vs. Fd curves for k equal to zero-drive state reciprocal 
lattice vectors (Fig. 0.a). This increase in order arise be- 
cause as Fd deacreses the vortices spend more of their 
travel time near the zero-drive state equilibrium posi- 
tions. There is no long range time order as shown by the 
absence of sharp peaks in the frequency dependence of 
C?rM andC^H ((Fig.|.b) 

iv) Plastic flow: low T, F d < F c , B = 1.25£ . For F d 
just below F c some spatial order remains. For x = and 
x = a c d, P( r ) (Fig.Q.b) consists of isolated spots, indicat- 
ing that the double vortex chain-structure of the infinite- 
drive state transforms to a single chain one with, roughly, 
nearest neighbor chains displaced relative to each other 
by half intra chain spacing, and that the vortices oscil- 
late around the corresponding equilibrium positions. For 
larger x the spots become interconnected by stripes in the 
direction of motion, indicating that inter chain relative 
motion at these separations, and thus plastic flow, is tak- 
ing place. The behavior of AVy m /Vy m across the transi- 
tion (Fig. ^.b) shows little change, suggesting that plastic 
flow is weak. The frequency dependence of C^ m (o;) and 
Off (u) is found to remain sharply peaked. These results 
are interpreted as indicating that vortices move together 
with the chain, oscillating around their equilibrium po- 
sition within it, but relative motion between chains is 
taking place. As Fd decreases further plastic flow contin- 
ues with increasing AVy m /Vy m , but the vortices remain 
ordered in single chains and are unpinned for Fd > 
(Fig. ^|.b). It is expected that vortex reordering and pin- 
ning would take place as Fd decreases. The failure to see 
this in the present simulations is interpreted as evidence 
that trapping in a metastable state is taking place. 
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FIG. 4. P(r) for B = 1.2550 in 512 x 512 space lattice: 
a) elastic flow phase at Fd — 1.3F P , b) plastic flow phase at 
F d = 1.2F P . 

v)T- dependence. As T increases F C (T) also increases, 
because thermal fluctuations make the infinite-drive state 
VL softer. Thermal fluctuations also lead to smaller val- 
ues of 5(k). This makes the dynamical transitions de- 



scribed above more difficult to observe, particularly for 
B = 1.25^0. For B = 2B^ the F C (T) curve is shown in 
Fig. |a. 

For other values of B > B^, and drive in the (0,1) or 
(l,0)-directions, the possible infinite-drive states are ei- 
ther a nearly triangular VL, or a VL with n vortex-chains 
trapped in each groove, so that a similar dynamical phase 
diagram is anticipated. 

The experimental verification of the new dynamical 
phases and of the transition predicted in this paper may 
be possible in clean superconducting films with sufficient 
thermal fluctuations to allow relaxation to the infinite- 
drive state to take place at large drives. Signatures of 
these phases are present in the V-I curve and in the noise 
power spectrum. Visualization of the spatial order may 
also be possible using decoration or fast STM techniques. 
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